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A simulation of quenched QCD with the overlap Dirac operator has been completed using 100 
Wilson gauge configurations at /3 = 6 on an 18 3 x 64 lattice and at /3 = 5.85 on a 14 3 x 48 lattice, 
both in Landau gauge. We present results for light meson and baryon masses, meson final state 
"wave functions," and other observables, as well as some details on the numerical techniques that 
were used. Our results indicate that scaling violations, if any, are small. We also present an analysis 
of diquark correlations using the quark propagators generated in our simulation. 
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I. INTRODUCTION 



Preserving chiral symmetry presented a notorious diffi- 
culty for the early formulations of lattice Quantum Chro- 
modynamics (QCD). The problem stemmed from the fact 
that an ultralocal discretization of the Dirac equation 
must either abandon full chiral symmetry or introduce 
extra fermionic degrees of freedom (doublers), with as- 
sociated no-go theorems for lattice fermions [1, 2]. It 
was only about a decade ago that seminal work by Ka- 
plan, Shamir, Neuberger, Narayanan, and others pro- 
vided a way out of the impasse, through the introduc- 
tion of domain- wall [3, 4] and overlap [5-7] formulations 
of lattice fermions. These two closely related discretiza- 
tions of the Dirac equation avoid the no-go theorems by 
forfeiting ultralocality, while retaining locality, and sat- 
isfy an identity, the Ginsparg- Wilson relation [8], which 
allows one to extend to the lattice the chiral symmetry 
of continuum QCD [9]. Unfortunately, the satisfactory 
properties of the new lattice formulations come with a 
heavy computational cost, since they entail either exten- 
sion of the lattice in a fifth dimension or very demanding 
large matrix manipulations. It is therefore important to 
subject these novel discretizations to the test of lattice 
QCD simulations on systems of realistically large size, in 
order to explore the adequacy of the necessary numerical 
techniques and to validate the good properties that are 
expected to follow from the preservation of chiral sym- 
metry. In this paper, we present the results of such an 
investigation, where we used the overlap Dirac operator 
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to simulate quenched QCD on two lattices, of size 18 3 x 64 
and 14 3 x 48. 

One early simulation of quenched QCD with overlap 
quarks was performed in [10, 11]. Pioneering investi- 
gations of QCD with the overlap fermion discretization 
were also presented in [12-23]. In the work of [10, 11], 
quenched QCD was simulated with the Wilson gauge ac- 
tion at j3 = 6 on a lattice of size 16 3 x 32. An expan- 
sion into fractions, after projection of a small number of 
low lying eigenvectors, was used to calculate the overlap 
Dirac operator. It was shown that the available com- 
putational methods could produce accurate results for 
the quark propagators in a reasonable amount of time 
with moderately powerful computer resources (capable 
of 10 to 20 Gflops sustained). Satisfactory results were 
obtained on the pseudoscalar spectrum, strange quark 
mass, quark condensate, renormalization constants, and 
a few other observables, and the good chiral behavior of 
the theory was verified. It also became apparent, how- 
ever, that the extent of the lattice in the temporal direc- 
tion was too small for a meaningful calculation of other 
parameters of the hadron spectrum, such as vector me- 
son and baryon masses, as well as for a reliable estimate 
of important matrix elements. We decided therefore to 
extend the scope of the investigation by considering a 
bigger lattice, with a time extent twice as large. The 
spatial extent was also slightly increased, with the final 
choice of lattice size, 18 3 x 64, motivated by a careful as- 
sessment of the computer resources we could rely upon. 
After completing the calculation of the quark propaga- 
tors on the 18 3 x 64 lattice, we also simulated a system of 
size 14 3 x 48, with correspondingly coarser lattice spac- 
ing, in order to check scaling. 

For the calculations described in this paper we simu- 
lated quenched QCD with Wilson gauge action at (2 = 6 
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on the 18 3 x 64 lattice and (3 = 5.85 on the 14 3 x 48 
lattice. The corresponding values of the lattice spac- 
ing are a -1 = 2.12 GcV and a -1 = 1.61 GcV, on the 
basis of the Sommer scale defined by r^Ffro) = 1.65, 
ro = 0.5 fm [24, 25]. (The determination of a -1 from the 
Sommer scale in quenched QCD is quite accurate and 
we can consider the corresponding statistical error neg- 
ligible with respect to the other statistical errors in this 
work.) Thus, the two lattices are of approximately the 
same physical volume. For both lattices, we generated 
100 gauge configurations using the multihit Metropolis 
algorithm, with 11,000 lattice upgrades for equilibration 
and 10,000 lattice upgrades between subsequent config- 
urations. For each configuration, we performed a gauge 
transformation to Landau gauge. We then calculated the 
quark propagators with source at the origin and all 12 
source color-spin combinations for bare quark mass 

am e {0.03, 0.04, 0.06, 0.08, 0.1, 0.25, 0.5, 0.75} (1) 

on the finer 18 3 x 64 lattice, and 

am G {0.03, 0.04, 0.053, 0.08, 0.106, 0.132, 0.33, 0.66, 0.99} 

(2) 

on the coarser 14 3 x 48 lattice. The calculation of the 
quark propagators was done with a conjugate-gradient 
multimass solver. Technical details of our calculation, as 
well as relevant formulae, are presented in Appendix A. 
We should mention here, however, some computational 
considerations which informed the choices we had to 
make for our investigation. 

The limitations of the quenched approximation are 
well-known, and forefront calculations today tend to in- 
corporate dynamical quarks. A simulation with dynam- 
ical overlap quarks on a lattice similar to the one we 
studied would have required, however, computational re- 
sources at least two orders of magnitude larger than we 
had available. Moreover, an efficient calculation of quark 
propagators with the overlap Dirac operator is a neces- 
sary prerequisite for overlap simulations with dynami- 
cal quarks. We concluded that a quenched calculation 
would be the best option at present, in that it would al- 
low us to consider a reasonably large lattice and separate 
the computational problems associated with the calcula- 
tion of the overlap operator from those of the dynamical 
fcrmion feedback. An alternative would have been to 
perform mixed action calculations, i.e. to use overlap va- 
lence quarks in the background of gauge configurations 
generated with some other type of dynamical quarks, for 
example those made available by the MILC collaboration 
(see http://qcd.nersc.gov). While that might have been 
computationally feasible and may be a good strategy for 
future calculations, we decided not to proceed in that 
manner at this stage in order to avoid the entanglement 
of two different sets of computational effects as well as for 
a more technical reason, on which we shall now elaborate. 

The overlap Dirac operator for a massless quark is 



given by 

where H(p) stands for the Hermitian Wilson-Dirac oper- 
ator with mass —p/ 'a, namely 

H(p) = y 5 D w (p) 

= f EMV,-Vt)+aVtV,]- ^ (4) 

with the forward lattice covariant derivative implicitly 
defined by 

V^{x) = - [U^(x)^(x + ap) - ip(x)] . (5) 

In terms of D, the overlap Dirac operator for a quark of 
mass m is then given by [1 — (am) / (2p)]D + m. (See Ap- 
pendix A for more details on the overlap Dirac operator.) 

In the calculation of the quark propagators, which 
is based on a conjugate-gradient iterative procedure, at 
each step one must apply the overlap operator to the cur- 
rent iterate. Equation (3) shows that this requires evalu- 
ating the action of (7? 2 ) -1 / 2 on a quark field, for which in 
turn one must use some suitable approximation to the in- 
verse square root. Implementing such an approximation, 
with the required degree of numerical accuracy, becomes 
more and more demanding the larger the condition num- 
ber of H 2 , i.e. the larger the value of the ratio between its 
largest and smallest eigenvalues. Since the largest eigen- 
value of H 2 is bounded, in practice the condition number 
of the matrix depends on its lowest eigenvalue, or, cquiva- 
lcntly, on the gap in the eigenvalues of H around 0. This 
gap in turn depends on the parameter p (which must 
be carefully chosen) as well as, loosely speaking, on the 
"amount of disorder" of the gauge configuration. This be- 
came quite apparent in our study, where the calculation 
of the quark propagators on the smaller 14 3 x 48 lattice 
turned out to be computationally more challenging than 
the calculation on the 18 3 x 64 lattice, because of the 
increased "roughness" of the background gauge configu- 
rations at the smaller value of (3. A similar dependence 
is encountered in the domain-wall formulation, where the 
so called "residual mass," which measures the deviation 
from chiral symmetry induced by the truncation in the 
fifth dimension, is seen to depend on the smoothness of 
the gauge configuration [26] . Given the fact that the gap 
in the eigenvalues of H can be substantially affected by 
the choice of action and by the presence or absence of dy- 
namical fcrmions, we thought that it would be useful to 
establish a benchmark by performing a quenched simula- 
tion with the most traditional action, namely the Wilson 
gauge action. We thought that demonstrating that the 
overlap formulation is amenable to precise calculations on 
large lattices in this context would be an important step 
toward the use of more elaborate gauge actions and/or 
mixed action calculations. 
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Another choice we had to make concerned the type of 
quark sources to use. Calculations of the lowest masses in 
the hadron spectrum are made more precise by the use 
of extended sources. However, we also wanted to take 
advantage of our quark propagators for the study of non- 
perturbative renormalization and the evaluation of se- 
lected matrix elements, both of which required point-like 
sources. Since the high cost of evaluating overlap quark 
propagators made it impossible for us to use both source 
terms, we had to adopt point-like sources. Of course, 
this did not prevent us from using extended sink opera- 
tors, and indeed some of the results we present here have 
been obtained with extended sinks. Finally the non- 
perturbative renormalization techniques required that 
the gauge field be brought to Landau gauge prior to 
evaluating the propagators. We have implemented this 
gauge fixing for all our gauge configurations (of course, 
after they have been generated), an additional advantage 
being that this allowed us to calculate quark-antiquark 
correlation functions in the final state as well as diquark 
propagators. 

Novel results with respect to the calculations pre- 
sented in [10, 11] have been obtained for the vector 
meson spectrum and decay constants, for the baryon 
spectrum, for correlations among quarks and antiquarks 
in the final state (which we will refer to as final state 
"wave functions"), scaling, quark and diquark propaga- 
tors, non-perturbative renormalization, and meson ma- 
trix elements of selected four-quark operators. In this 
paper, we present our results for the meson spectrum 
and decay constants, meson final state wave functions, 
the baryon spectrum, scaling from (3 = 5.85 to /3 = 6, 
and quark and diquark propagators. We will present our 
results for renormalization constants and matrix elements 
in a companion paper. We hope to present a more de- 
tailed analysis of heavy quark states, as well as a study 
of diquark correlations within baryons, in future publi- 
cations. 

In Section II, we present our results for light meson 
observables on the finer 18 3 x 64 lattice. We follow in 
Section III with our results for the coarser 14 3 x 48 lat- 
tice, as well as a scaling comparison between the two 
lattices. As discussed in Appendix A, we chose p = 1.4 
at = 6.0 and p = 1.6 at f3 = 5.85 [13, 27] so as to 
optimize the locality properties of our overlap operator. 
That choice induces a lattice spacing dependence which 
would have to be parameterized in order to allow for a 
continuum extrapolation. Since our goal is only an esti- 
mate of possible discretization errors, we do not pursue 
such an approach here. In Section IV, we present our 
results for the baryon spectrum on the 18 3 x 64 lattice 
as well as scaling between the two lattices. Finally, in 
Section V, we present our analysis of diquark correlation 
functions and diquark spectra. In the Appendices, we 
give background information on the overlap Dirac oper- 
ator, present details of our simulation techniques, and 
present data tables for meson and baryon spectra and 
other observables. 
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FIG. 1: Time dependence of a PP correlator, for quark 
masses ami = 0.08, «ni2 = 0.10. 



II. LIGHT MESON OBSERVABLES ON THE 
18 3 x 64 LATTICE 

A. Correlators 

We evaluated meson correlators for point source, point 
sink operators, as well as point source, extended sink op- 
erators. As explained in the introduction, computational 
limitations did not allow us to calculate correlators with 
extended sources. 

We first consider correlators with point sink operators. 
We only consider connected diagrams, so without loss of 
generality we take quark and antiquark of different flavors 
fi and /2 . The zero-momentum meson correlator is given 
by 

Gab® = (^Tr[^ 2 (0;x,/)r A7 5(5 /l (0;^/)) t 75r B ]) , 

X 

(6) 

where S(0\x,t) is the Euclidean quark propagator and 
Ta and Tb are the Dirac 7-matrix combinations associ- 
ated with the meson states A and B. An example of such 
a meson correlator is shown in Fig. 1. 

To obtain ground state meson masses, the zero- 
momentum meson correlators, which are even under 
time-reversal, were fit in an appropriate fitting window 
to the usual functional form 



G(t) 



Z 
M 



e- MT / 2 cosh 



(7) 



where M is the meson mass and T is the extent of the 
lattice in time, and where we will refer to Z as the cor- 
relator matrix element. 

For mesons with quarks of non-degenerate mass, we 
combined the results obtained by an interchange of the 
two quark masses. In order to increase statistics, we 
folded the meson correlators over the midpoint T/2 and 
fitted to data below that point. Concerning statistical 
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FIG. 2: Effective mass plateau for different pseudo-scalar cor- 
relators at am = 0.08. 



FIG. 3: Scanning for optimum t m in for PP mass, with ami 
0.08, airi2 — 0.10, and t max = 32. 



errors, the data that enter into the fit (the meson corre- 
lators at different values of time) are highly correlated, 
but the amount of data is not sufficient for an evaluation 
of cross-correlations that would be precise enough for in- 
corporation into the fitting procedure. Thus, we instead 
used the bootstrap technique. We found that the error 
estimates became stable when the number of bootstrap 
samples jib reached a value of approximately 200 and 
used tib = 300 in our calculations. We used the boot- 
strap method for the estimate of the statistical errors in 
most of the results presented in this paper. 

Typically, effective mass plots are used to determine 
the best fitting window {t m in,t max } for correlators. Ex- 
amples of such effective mass plots are shown in Fig. 2. 
Using a cosh fit for the meson correlators allowed a con- 
sistent t max of T/2. In order to pinpoint the best value 
for t m in , we considered scans over different values of t m i n 
for a fixed t max . An example of such data is shown 
in Fig. 3. The value of t m i n chosen was the smallest 
value (consistent with the errors) before the clear effect 
of higher states caused the mass prediction to rise. Con- 
sideration of the x 2 value from the fit also was used to 
confirm the choice. If possible, a single value of t m i„ was 
chosen for all quark mass combinations. 



B. Meson spectra 

The effective mass data and fitting window scans in 
Figs. 2 and 3 suggested use of a fit range 12 < t/a < 32 
in order to extract the pscudoscalar mass. For the vector 
mass, the fitting range was 8 < t/a < 32. We illus- 
trate in Fig. 4 our results for the pseudoscalar spectrum 
for all possible input quark mass combinations. In the 
quenched approximation, the correlator Gpp{t) receives 
contributions proportional to 1/m 2 and 1/m from chi- 
ral zero modes that are not suppressed by the fermionic 
determinant. These unsuppressed contributions, which 
should vanish in the infinite volume limit, could be siz- 



1.5 



8 



S 



0.5- 



O PP point sink 

□ PP ext. sink 

V PP-SS point sink 



0.2 0.4 0.6 

am = a(m 1 +m 2 )/2 



0.8 



FIG. 4: Pseudoscalar meson spectrum for point and extended 
sinks. 



able at finite volume [13, 26]. A method of handling 
these quenching artifacts is to consider the difference of 
pseudoscalar and scalar meson correlators 



G 



PP-SS 



(t) = Gpp(t)-G ss (t) 



(8) 



since, by chirality, the quenching artifacts cancel in the 
difference [26]. The results for the pseudoscalar masses 
obtained with this correlator are also included in Fig. 4. 
With our large lattice we do not see any significant differ- 
ence between the results obtained with PP and PP— SS 
correlators. Thus, in the remainder of this paper, for 
what concerns the pscudoscalar masses and matrix ele- 
ments, we will only consider the PP correlators, except 
where otherwise noted. For economy of figures, we also 
show in Fig. 4 the results for correlators with extended 
sink operators, which will be discussed in section II D. 

The chiral behavior of the pseudoscalar spectrum with 
degenerate quarks is illustrated in Fig. 5 for the PP and 
PP — SS correlators. The lightest data point corresponds 
roughly to a kaon composed of degenerate quarks of mass 
am = 0.03 ~ am s /2, where am — a{m\ + m-i)/'! is the 
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FIG. 5: Chiral behavior of meson masses. The solid line 
corresponds to the best fit to Eq. (9) in the interval am < 0.1. 



average of the quark and antiquark masses ami and ami 
in the meson and am s is the bare strange quark mass, 
all in lattice units. We neglect for the moment the chiral 
logarithms discussed in Section II E and perform a fit to 
the linear form 



FIG. 6: Vector meson spectrum for point and extended sinks. 



{aMpf = A + B(am) 



(9) 



where aMp is the pseudoscalar mass. If we fit all data 
with am < 0.1 we obtain (see Fig. 5) 



A = 0.0058(15) , 
for the PP correlator and 
A = 0.0059(16) , 



B = 1.376(15) 



B = 1.380(17) 



(10) 



(11) 



for the PP — SS correlator. (The corresponding val- 
ues in [10] were A = 0.006(4), B = 1.39(3) and A = 
-0.0005(68), B = 1.43(7) for the PP and PP-SS corre- 
lators, respectively.) These results exhibit the good chiral 
behavior of the overlap formulation. The non-zero value 
of the intercept A in Eqs. (10) and (11), albeit small, 
is however statistically significant. Since our results rule 
out the possibility that this is due to zero modes, the 
small deviation from chiral behavior can originate cither 
from finite volume effects or chiral logarithms. In a later 
section we will argue against finite volume effects and 
show that it is indeed compatible with chiral logs. 

Using a larger lattice than [10] allowed for observation 
of vector meson states, as shown in Fig. 6, with aMy 
the vector meson mass. The observation of vector meson 
states was enhanced by use of extended sink operators, 
as discussed in Section II D. Our results for the meson 
spectrum are reproduced in the tables presented in Ap- 
pendix B. 



C. Axial Ward identity and Za 

Exact chiral symmetry implies a conserved axial cur- 
rent, and the associated axial Ward identity (AWI) prc- 
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FIG. 7: AWI ratio as a function of time for all degenerate 
quark mass combinations. 



diets a constant value for the ratio 



Pit) 



G PP (t) 



(12) 



The conserved axial current is a local, but not ultralocal, 
operator. The ultralocal axial current 



Aq = -0i(z)7o75 



;i- w-D)^ 

2p 



(x) 



(13) 



differs from the exactly conserved axial current by a fi- 
nite renormalization factor Za and possible corrections 
0(a 2 ). We calculated the correlator in Eq. (12) with the 
current of Eq. (13) and using the lattice central differ- 
ence for Vo, corrected so as to take into account the sinh 
behavior of the correlator. Figure 7 shows the observa- 
tion of plateaus for all of our quark masses in the range 
8 < t/a < 56. 

The fit shown in Fig. 8 to 



ap = A + 2(am)/ZA + C(am) 



(14) 
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FIG. 9: PP extended sink correlators G(r,t) for various t/a 
and ami = am.2 = 0.03. 



gives 

A = 0.00002(10) , Z A = 1.5555(47) , C = 0.273(32) . 

(15) 

(In [10], the results were A = -0.00002(7), Z A = 
1.555(4) and C = 0.277(12).) The fact that the value 
of A is consistent with zero is an excellent indication of 
the good chiral behavior of overlap fermions (compared 
to the residual mass found in domain-wall fermion cal- 
culations). Also, C is rather small, a possible indication 
that discretization errors might be smaller than expected 
on the basis of purely dimensional arguments. 



D. Extended sink operators 

For some correlators, e.g. the vector VV correlators, 
the signal for the ground state is small when point sources 
and point sinks are used. In particular, the use of point 
sources and sinks causes coupling to excited states that 
do not decay until relatively large values of t m i n . The sig- 
nal can be improved by using non-local extended opera- 
tors instead of local point operators in the representation 
of the meson [28]. Since, as mentioned in Section I, our 
quark propagators were calculated using point sources, 
we were limited only to the case of extended sinks. Here 
we elaborate on our use of extended sink operators. 

We consider the correlator 

G(r,t) = (Tr[S^(0;x,t)T Al5 (S fl (O;y,0 l5 T B ] 
x S(\x-y\-r)) , (16) 



■r.y 



where the ^-function is approximated on the lattice with 
a Kronccker 5, taking into account the multiplicity of the 
sites. The quantity r is the separation of the quark and 
antiquark at the sink in lattice units. No gauge trans- 
port factor is needed to make the non-local correlator in 
Eq. (16) well defined because we calculated the quark 
propagators in Landau gauge. 
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FIG. 10: Comparison of t m in scans for VV point and extended 
sinks with ami = ami = 0.03. 



We calculated G(r, t) for PP and VV correlators. 
A fast Fourier transform (FFT) was used in order to 
speed up the calculation; we first Fourier transformed 
the quark propagators and then used the convolution the- 
orem. That reduced the double summation over spatial 
lattice sites to a single sum, decreasing the computational 
time by almost three orders of magnitude [29]. Figure 9 
illustrates the average value of the G(r, t) correlators for 
pseudoscalar mesons with quark mass am = 0.03, each 
normalized to unity at r = 0. A clear ground state "wave 
function" is apparent after about t/a = 8. 

The observation that the wave function settles into 
a definite profile representing the contribution of the 
ground state allows us to use this ground state wave 
function to construct extended sink correlators. For both 
the PP and VV correlators, we used the corresponding 
function tp(r) = G(r, 8a) (of course, with different <p(r) 
for the pseudoscalar and vector states) to define an ex- 
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FIG. 11: Vector meson spectrum for am < 0.1. 



FIG. 12: Evidence for quenched chiral logs. 



tended sink correlator 

G ext (t)=5>(r)G(r,i) 



(17) 



from which we then extracted the meson mass. 

The use of extended sinks was most valuable in the cal- 
culation of the vector meson spectrum, due to the result- 
ing increase in the size of the fitting window (see Fig. 10). 
The point sink and extended sink vector meson spectra 
are compared in Fig. 11 for low quark mass, am < 0.1. 
A linear fit of the quark-mass dependence of the vector 
meson mass obtained from the extended sink data gives 



aM v = 0.409(15) + 1.10(13) (am) . 



(18) 



We note that the chiral limit value of 0.409(15) is larger 
than the value 0.366 obtained using the Sommer scale 
value for the lattice spacing and the experimental p mass. 
As we will discuss in a later section, the discrepancy 
gives an indication of the systematic error induced by 
the quenched approximation. 



E. Quenched chiral logarithms 

A linear fit to the pseudoscalar masses for quark mass 
am < 0.1 (see Fig. 5) produces a line with an inter- 
cept which is very small, but statistically different from 
zero, indicating a deviation from the chiral behavior 
(aMp) 2 (x am. This deviation is put in better evidence 
by considering the ratio (aMp) 2 / am, shown in Fig. 12, 
which exhibits a sharp rise at low am. 

Having already ruled out that the non-zero intercept 
may be due to the effect of fermion zero modes (the PP 
and PP — SS correlators produce statistically indistin- 
guishable results for the intercept), the small deviation 
from linear behavior could be due to finite volume ef- 
fects or to chiral logs. Insofar as finite volume effects are 
concerned, not having the resources needed to repeat the 
calculation on a larger lattice, we can only observe that 



the Compton wavelength for our lightest pseudoscalar, 
Alp 1 ~ 4.5a, is much smaller than the size of our lattice, 
L = 18a (AIpL = 4). In the rest of this section, we com- 
pare our results for small quark masses with the predic- 
tions from quenched chiral perturbation theory [30, 31]. 

We fit the degenerate quark mass results for the pseu- 
doscalar masses to the expression [30] 



(aM F 



A(am) 



1/(1+5) 



B(am) 



(19) 



where the leading quenched logarithms proportional to 
S = 2m 2 /(N c (4irf x ) 2 ) have been resummed into a power 
behavior and where the term proportional to B param- 
eterizes possible higher-order corrections in the mass ex- 
pansion. Here mo is the singlet contribution to the rf 
mass, N c = 3 the number of colors, and f x the value of 
the pion decay constant in the chiral limit. The fit is 
shown in Fig. 12, and the results for the parameters are 

A = 0.680(68), 5 = 2.98(31), (5 = 0.29(5), (20) 



of S 



presented elsewhere in the 

2 



consistent with values 

literature [26, 32-37]. Note that a fit to (aMp)" m 
Eq. (19) containing an additional constant term C does 
not change the central values of A, B, or 6 and produces 
a value for C that is very small and consistent with zero. 

For the non-degenerate quark case, a value for the 
quenched chiral log parameter 5 was obtained via a fit 
to the expression y = 1 + 6x [31, 34], where 



2m? 



M 2 

mi + ?7i2 M 2 X mi + m-i M| 2 



2mi M 2 , 



mi + m 2 

mi — 777,2 



In 



777 2 
7771 



(21) 



(22) 



and Mij is the mass of a pseudoscalar meson composed 
of quarks with masses and mj. In the derivation 
of this functional form, the contributions of the kinetic 
term of the singlet Lagrangian with coupling a, which is 
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FIG. 13: Calculating the quenched chiral log parameter 5 
with non-degenerate quark masses. 



sublcading in a 1/N C counting, were neglected, as they 
were in Eq. (19). However, here the quenched chiral log 
proportional to 6 is considered a correction and dealt with 
at linear order. In addition, all higher order terms in the 
chiral expansion are neglected. The value of S from our 
non-degenerate data is 



8 = 0.18(8) , 



(23) 



which is consistent, within the large statistical errors, 
with the result in Eq. (20). 

As shown in Section II B, these chiral logarithms in- 
duce only a small deviation from linear behavior in the 
relation of (aMp) 2 vs. (am) at the values of light quark 
mass used in our simulation. Moreover, our statisti- 
cal errors are still large and our lightest pscudoscalar 
meson has a mass around that of the kaon, where the 
quenched theory is tuned to reproduce the unqucnched 
theory. Therefore, we take the slope parameter B, ob- 
tained from the linear fit of Eqs. (9) and (10), to be our 
estimate of the coefficient of the leading term in the chiral 
expansion of {aMp) 2 as a function of quark mass am. 



F. Decay constants and determination of the 
lattice spacing 

We extract the decay constant fp from the relation 
2to|(0|P|7t)| = 2mVz = f P Mp , (24) 



where the correlator matrix clement Z is defined in 
Eq. (7). Our results for the pscudoscalar decay con- 
stant are reproduced in Fig. 14, where we plot afp as 
a function of (aMp) 2 , together with a line representing 
the results of a linear fit, as predicted by NLO quenched 
chiral perturbation theory [30, 31]. Following [38], we 
also determine the lattice spacing by the method of lat- 
tice physical planes. In particular, we are looking for 
the point in our lattice parameter space, at which the 
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FIG. 14: Results for fp and determination of a from the 
method of physical planes. The dashed line corresponds to a 
linear fit of the data for afp in the interval am < 0.1; the solid 
curve represents the parabola (afp)/(aMp) — (/jr/Mjf)exp. 



dimcnsionless ratio of the pscudoscalar mass and decay 
constant is equal to the experimentally determined ratio 
of the Kaon mass to decay constant. The parabola in 
Fig. 14 shows the line, along which 



afp 
aMp 



(JjL 



0.323 , 



(25) 



where (/k/Mk)cxp stands for the experimental value 
of this ratio and we used as experimental data Mk = 
0.495 GcV and f K = 0.16 GeV [39]. The intersection of 
the two lines gives 



aM K = 0.226(6), af K = 0.074(2) 



(26) 



from which, using the input Mk = 0.495 GcV, one gets 

a^ 1 = 2.19(6) GeV . (27) 

Using the Sommer scale value of the lattice spacing 
and the experimental value of the pion mass, M w = 
0.135 GeV [39], yields 



and thus 



af„ = 0.065(2) 



f -f- = 1.13(4) , 



(28) 



(29) 



which is below the experimental value of 1.22 but is com- 
patible with other quenched calculations [40-44]. 

The method of lattice physical planes has the advan- 
tage of using data in a region of quark masses accessible 
to the lattice calculation, avoiding the need to perform a 
chiral extrapolation to low quark mass. 

The value for a -1 derived with the method of lattice 
physical planes should be contrasted with the one derived 
from the Sommer scale, namely a -1 = 2.12 GcV at (3 = 



9 



6. A further independent determination of the lattice 
spacing can be obtained by extrapolating our results for 
the vector meson spectrum to the chiral limit. This gives 



0.2 r 



aM p 
infer 



0.409(15) (see Eq. (18)), from which one would 



a- 1 = 1.90(4) GeV 



(30) 



The discrepancy among the three values of a obtained 
above, namely a -1 = 1.90(7) GeV from the p mass, 
a -1 = 2.12 GeV from the Sommer scale, and a -1 = 
2.19(6) GeV from the method of physical planes, is sub- 
stantially larger than what could be due to statistical er- 
rors alone and should be attributed for the most part to 
the quenched approximation. The other sources of error, 
namely those due to the finite lattice spacing, extrapo- 
lation to small quark masses, and finite volume effects, 
are substantially smaller, as can be inferred from the data 
presented in this paper or, insofar as finite volume effects 
are concerned, argued from the size of the lattice. Tak- 
ing the maximum variance in the three numbers above, 
Aa _1 = 0.29 GeV, as an indication of the systematic er- 
rors in the quenched approximation, the corresponding 
relative error is Aa _1 /a _1 ~ 14%. In the rest of this pa- 
per, whenever we quote quantities in physical units rather 
than lattice units, we will use the value of a -1 from the 
Sommer scale for the conversion and systematic errors on 
such quantities, when given, will include an estimate of 
the scale setting ambiguity from the method of physical 
planes. 

Equation (26) together with the linear fit of Eqs. (9) 
and (10), which was justified at the end of Section HE, 
gives us 



a(m s 



= 0.0661(44) 



(31) 



where m s stands for the bare mass of the s-quark and rh 
for the average bare mass of the light u and d quarks. We 
note in passing that using the value of Eq. (31), together 
with our results for the vector meson masses (from the 
extended sink correlators), give 



M, 



K ' 



= 1.09(5) 



(32) 



which is compatible within errors with the experimen- 
tal value of 1.15 as well as with other quenched calcu- 
lations [40]. Use of the Sommer scale value of the lat- 
tice spacing, together with our data for the pseudoscalar 
spectrum and the experimental value for the kaon mass, 
produces a slightly larger value for a(m s + rh) than the 
one obtained from the method of physical planes, namely 



a(m s +m) = 0.0709(17) 



(33) 



We will use this value in the rest of the paper. The vari- 
ation that the change in a(m s + rh) induces in the values 
already quoted for fa/ fir and Mk* /M p are minimal and 
well below the statistical errors. 
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FIG. 15: Vector decay constant vs. bare quark mass. 



We calculated the vector decay constant using 



ah 



aZ\, 



My 



(34) 



where Z is the matrix element appearing in the cosh fit 
of the point source, point sink vector meson correlator 
and Zy is the renormalization constant of the ultralocal 
vector current. For My, we used the extended sink vec- 
tor masses. Like the axial current, the conserved vector 
current in the overlap discretization is a local, but not ul- 
tralocal, operator, and the ultralocal current used in the 
vector meson correlator must be renormalized. The chiral 
symmetry properties of the overlap formulation guaran- 
tee, however, that Zy = Za, and so we can use the result 
obtained in Sect. II C. 

Unfortunately, as described in Sect. II D, the the cou- 
pling of the vector meson ground state to the point sink 
is rather small. Therefore it is difficult to find a valid 
plateau range especially for small quark masses. In fact, 
as can be seen in Fig. 11, only at around am > 0.08 
do the masses extracted from point and extended sink 
propagators coincide. We therefore also expect the Z 
extracted from smaller masses to be heavily affected by 
excited state contributions and choose to ignore them. 
Consequentlly, we only use the points with quark mass 
am = 0.08, 0.10, 0.25 in the fit of the afv data. Our re- 
sults for afv are shown in Fig. 15. A linear fit of the 
data, along with the quark mass of Eq. (33), yields 

af p = 0.125(5), af K , = 0.128(5), f -f- = 1.03(6) . 

Jp 

(35) 

The ratio fx* / f P agrees well with the experimental value 
1.03(4) [40]. 

Using the Sommer scale value for the lattice spacing, 
the above results give in turn 



f p = 265(11) MeV, f K , = 272(10) McV 



(36) 
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G. Quark masses and chiral condensate 

Our result for the bare quark masses, a(m s +m), can be 
converted into a corresponding result for the renormal- 
ized quark masses. The bare quark mass m(a) is related 
to the renormalized quark mass m([i) by 



m{ii) = lim Z m (afi)m(a) 



(37) 



The mass renormalization constant Z m is in turn related 
to the renormalization constant Zs for the non-singlet 
scalar density by Z m (a\x) = 1/Zs(afj,). We calculated 
Zs in the RI-MOM scheme starting from the identity 



Zjp(a^) = lim Z A 



T A (p, m) 



m^o Ts{p,m) 



(38) 



p z — fi z 



where Ta(p, m) and T$(p, m) are suitably defined 
Green's functions for the axial current and the scalar den- 
sity in Landau gauge and Za is the renormalization con- 
stant for the axial current calculated in Sect. II C. Details 
of the procedure will be presented in a separate publica- 
tion. The Green's functions in Eq. (38) have been calcu- 
lated non-perturbativcly in a window of momenta which 
extends into the perturbative QCD domain, where con- 
tact can be made with perturbatively calculated renor- 
malization constants. We extracted our central value 
of Zs by performing a combined fit to am = 0.03 and 
am = 0.1 data in a momentum range p 2 = 3 — 14GeV 2 
using four-loop running [45] and additional (ap) 2 terms 
to account for discretization effects as well as 1/p 2 and 
1/p 4 terms to account for mass effects and other possible 
subleading terms in the operator product expansion of 
the relevant correlation function. We obtain 



Z^(2GeV) = 1.25(2)(2) 



(39) 



where the first error is statistical and the second is an 
estimate of the systematic error obtained by varying the 
fit range and dropping additional terms when indicated. 
With this result, one can use the three-loop perturbative 
calculation of the ratio Z^/Z^ 1 [45] to calculate 

Zf S (2GeV) = 1.44(2)(3) . (40) 
Putting together Eqs. (33) and (40), we obtain 

(m s + m) m (2 GeV) = 105(3)(4) MeV (41) 

for the sum of strange and light quark masses. Using the 
value m s /m — 24.4(1.5) from chiral perturbation the- 
ory [46], we obtain 



MS 



(2 GeV) = 101 (3) (4) MeV 



(42) 



for the strange quark mass, which agrees well with the 
quenched lattice world average [41, 47, 48]. 

In the unquenched theory, the bare chiral condensate 
with valence overlap quarks is defined as 



X(a) = lim 



o N f 



^(0) 



2p 



D Up 



(0) 



(43) 



where m is the common mass given to the light quarks. 
It satisfies the integrated non-singlet chiral Ward identity 



1 

w f 



(0)Um^(P(x)P c (0)), 

' X 

(44) 

where P is the pseudoscalar density composed of two 
mass-degenerate quarks of different flavor and P c is the 
density obtained by interchanging the quark flavors. In- 
serting a complete set of states in (P(x)P c (0)} gives 



X(a) 



hm -^|(0|P|P)| 2 

ro— >0 Alp 



(45) 



where Mp is the mass of the pseudoscalar state \P). Us- 
ing 

2m\(0\P\ir)\=f P M P , (46) 

where fp is the pseudoscalar decay constant, yields the 
familiar Gell-Mann-Oakcs-Renner (GMOR) relation 



x{a) = - lim 



fpM 2 p 



m— >o 4m 



(47) 



Due to quenched chiral logarithms, x( a ) IS ill-defined 
in quenched QCD. However, as argued in Section HE, 
the slope parameter B, obtained from the linear fit of 
Eqs. (9) and (10), should be a reasonable estimate of 
lim m ^o Mp/m. Thus, we take 



1 



X (a) = --flBa- 1 



(48) 



as our determination of the physical quark condensate. 
Using our data for fp we get 



a 3 x(a) = -0.00144(10) 



X(a) = -0.0137(10) GeV 3 



(49) 



(50) 



Note that using the value f x = 0.123 GeV [49] would 
instead give X (a) = -0.0110(1) GeV 3 . 
Using $il>} m = Zf^xia), we finally get 

(^V> M ^(2GeV) = -0.0197(14)(20)GeV 3 

= -[270(6) (9) MeV] 3 . (51) 

Our results for m s and (ipip) are i n good agreement 
with the results presented in [10]. The values were 
m^(2GcV) = 102(6)(18)MeV and (ipip} m (2 GeV) = 
-[267(5)(15)MeV] 3 . They also agree with the determi- 
nations of the condensate using finite-size scaling tech- 
niques [13, 16], as well as with the recent continuum- 
limit calculation of [50] , all of which were obtained using 
overlap fermions. 
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FIG. 16: Direct determination of the chiral condensate. The 
curve corresponds to a quadratic fit of the data according to 
Eq. (54) in the interval am < 0.1. 



One could also attempt a determination of \ directly 
from a fit to the mass dependence of the quantity 



-a 3 x(m) = amJ2( P ( x ) pC ( )) 



(0) ) . (52) 



This is, however, made difficult by the very steep depen- 
dence of x(m) on m and also by possible infrared di- 
vergent contributions from zero modes. Contrary to the 
case of the determination of the pscudoscalar spectrum, 
where we found the effect of zero modes to be suppressed 
because of the large size of our lattice, zero modes are 
likely to contribute to the expression in Eq. (52), be- 
cause it involves the short distance behavior of the quark 
propagator. The contribution from zero modes can be 
eliminated by instead considering the subtracted expres- 
sion 

-a 3 X (m) = amJ2((P(x)P c m ~ (S(x)S c (0)}) , (53) 

X 

which has the same m — ► limit as am '}2 x {{P{x)P c (Q)) 
and where the contributions from chiral zero modes can- 
cel. Of course, in the quenched theory, that quantity 
diverges in the chiral limit due to quenched chiral loga- 
rithms. However, here again the effect of the quenched 
chiral logarithms appear to be small for the light quark 
masses reached in our simulation, and we assume that a 
polynomial extrapolation of our quenched results to the 
chiral limit gives a reliable estimate of the physical con- 
densate. Our results for — a 3 x(m) are shown in Fig. 16. 
A quadratic fit 



-a 3 x(m) 



-a 3 X 



B(am) + C(am) 2 



(54) 



gives the result 



-a 3 x = 0.00131(8) , B = 0.806(3), C = -0.14(1) 



(55) 
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FIG. 17: Pseudoscalar effective mass plateau for different 
channels at am — 0.106. 



The value we obtain in this way for a 3 x is compatible 
with the value in Eq. (49). The renormalized value is 



(V>V> Mb (2GcV) 



-0.0179(ll)(18)GcV 3 
-[262(5) (9) Me V] 3 . (56) 



III. MESON SCALING ANALYSIS 

Here we present our results for the coarser 14 3 x 48 
lattice, with (3 = 5.85, as well as comparisons between 
the two lattices. 



A. Meson spectra 

We extract meson masses from the correlators Gpp(t) 
and Gpp-ss(t)- The effective mass plateau is plotted in 
Fig. 17 for quark mass am = 0.106. The data suggest 
use of a symmetrized fit range 10 < t/a < 24 in order to 
extract the pseudoscalar mass. 

The extracted meson masses and matrix elements us- 
ing this fit range are reported in Tables VI and VII in Ap- 
pendix B. We neglect for the moment chiral logarithms 
and perform a fit to 



(aM P ) 2 =A + B(am) . 



(57) 



Fitting all data with am < 0.132, we obtain (see Fig. 18) 
^4 = 0.0045(15), B= 1.923(16) (58) 
in the PP channel and 



A = 0.0007(29) , B = 1.959(23) 



(59) 



in the PP — SS channel. In general, as for the 18 3 x 64 
lattice, the PP and PP — SS channel results were com- 
patible on the 14 3 x 48 lattice, so all further pscudoscalar 
data, unless otherwise specified, arc from the PP chan- 
nel. 
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FIG. 18: Chiral behavior of meson masses. The lines corre- 
spond to the best fit to Eq. (57) in the interval am < 0.132. 
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FIG. 20: The plateau of the AWI ratio ap(t) for different bare 
quark masses. 
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FIG. 19: Vector meson masses on the coarser lattice. The 
line corresponds to the best fit of the extended sink data to 
Eq. (60) in the interval am < 0.132. 



FIG. 21: The plateau value of ap vs. bare quark mass. The 
curve corresponds to the best fit to Eq. (14) in the interval 
am < 0.132. 



A fit of the quark-mass dependence of the vector meson 
mass obtained from the extended sink data, shown in 
Fig. 19, gives 

aM v = 0.605(29) + 0.84(22) (am) . (60) 

We note that the chiral limit value of 0.605(29) is larger 
than the value 0.483 obtained using the Sommer scale 
value of the lattice spacing and the experimental p mass. 
As for the finer lattice, we expect that this is a result of 
quenching errors. 



discretization effects, since that quark mass is large and 
comparable to the inverse lattice spacing. 

Taking the aforementioned plateau range and perform- 
ing a fit toEq. (14) with all data for am < 0.132 (Fig. 21), 
we obtain 

A = 0.00004(5) , Z A = 1.4434(18) , C = 0.381(8) . 

(61) 

Our results are compatible with chiral symmetry for all 
bare quark masses except am = 0.99 and the resulting 
Za is in agreement with [51, 52]. 



B. Axial Ward identity and Za 

Figure 20 shows the ratio p(t) = Gy A p(t)/Gpp(t) 
for all of our bare quark masses. We observe a nice 
plateau for all quark masses, except the largest bare 
quark mass am = 0.99, in a range 8 < t < 40. The 
lack of a plateau for am = 0.99 we believe is due to 



C. Quenched chiral logarithms 

A fit of the degenerate quark mass PP channel results 
to the expression (aM P ) 2 = Aiam) 1 ^ 1 ^ + B(am) 2 of 
Eq. (19) (sec Fig. 22) gives 

A = 1.20(11) , 5 = 2.66(41), 5 = 0.17(4). (62) 
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FIG. 22: Evidence for quenched chiral logs. 
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FIG. 23: Calculating the quenched chiral log parameter S. 



Considering the case of unequal quark masses, a fit to 
the expressions of Eqs. (21) and (22) (see Fig. 23) gives 
6 = 0.22(4). 

As discussed in Section HE for the results obtained 
on the finer lattice, we will take the slope parameter B, 
obtained from the linear fit of Eqs. (57) and (58), to be 
our estimate of the coefficient of the leading term in the 
chiral expansion of (aMp) 2 as a function of quark mass 
am. 



D. Decay constants and determination of the 
lattice spacing 



As for the 18 3 x 64 lattice, we plot in Fig. 24 
the lattice decay constant afp versus meson masses 
(aMp) 2 as obtained by the simulation. The continuous 
curve represents the physical ratio of these quantities, 
= 0.323. It turns out that the mesons at 



(fa/mic) ex 
our bare quark mass am 
to the physical kaons. 
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FIG. 24: Results for fp and determination of a from the 
method of physical planes. The dashed line corresponds 
to a linear fit of the data for afp in the interval am < 
0.132, the solid curve represents the parabola (afp) /(aMp) = 
(/k/MrtW 



The intersection of the two curves gives 

aMfc = 0.343(9) , af K = 0.109(2) , (63) 

from which, using the input Mk = 0.495 GeV, one gets 

a' 1 = 1.44(4) GcV. (64) 

Use of the Sommcr scale value of the lattice spacing 
and the experimental pion mass yields 



af w = 0.100(3) 



and thus 



fx 



1.09(4) , 



(65) 



(66) 



which agrees with Eq. (29). 

As for the 18 3 x 64 lattice, the value for a -1 derived 
with the method of lattice physical planes should be con- 
trasted with the one derived from the Sommer scale, 
namely a -1 = 1.61 GeV at p = 5.85. The lattice spacing 
obtained by extrapolating our results for the vector me- 
son spectrum to the chiral limit, aM p = 0.605(29) (see 
Eq. (60)), is 



= 1.28(6) GeV 



(67) 



As for the finer lattice, the discrepancy between the 
three values for the lattice spacing can be attributed 
to the quenched approximation, with a relative error of 
Aa-Va- 1 ~ 20%. 

Equation (63), together with the linear fit of Eqs. (57) 
and (58), which was argued for at the end of Section III C, 
gives 



a(m s +m) = 0.1176(59) 



(68) 
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FIG. 25: Vector decay constant vs. bare quark mass. 

As for the finer lattice, we note in passing that the value 
of Eq. (68), together with our results for meson masses, 
give 



M 



K ' 



Mo 



= 1.07(6) 



(69) 



Using the Sommer scale value of the lattice spacing and 
our pseudoscalar spectrum data yields a slightly smaller 
value for a(m s + to) of 



a(m s 



= 0.0943(15) 



(70) 



We will use this value for the rest of the scaling analysis. 

We again calculat ed the vector decay constant using 
afy = aZyy/Z/My. In the fit of the afv data, we used 
only the points with quark mass am = 0.106, 0.132, 0.33, 
since the point-point correlators give poor results for Z 
for lower quark masses. Our results for afv are shown 
in Fig. 25. A linear fit of the data, along with the quark 
mass of Eq. (70), yields 



af p = 0.175(13) , af K * = 0.179(12) , ^f- = 1.02(10) . 

fp 

(71) 

The ratio /#* / f p agrees well with both the value on 
the finer lattice, Eq. (35), and the experimental value 
1.03(4) [40]. 

Using the Sommer scale value for the lattice spacing, 
the above results give in turn 

f p = 281(22) McV , f K , = 287(19) McV . (72) 



E. Quark masses and chiral condensate 

As in Section II G, we calculate Z s in the RI-MOM 
scheme from 



Zf{ af x)= limZ. 



Ta(p, to) 



m— »o Ts(p,m) 



(73) 



Only data with am < 0.132 were used for this purpose. 
In order to obtain our final result, we used a combined 
fit to am = 0.04 and am = 0.132 in the range p 2 = 2 — 
8 GeV 2 for central values and obtained an estimate of the 
systematic error by varying the fit range and dropping 
additional terms when indicated. We fit Zs to the same 
functional form as for (3 = 6.0, yielding 



Z s u {2 GeV) = 1.29(3)(14) 



and 



/MS 



(2 GeV) = 1.48(3)(16) 



(74) 



(75) 



Putting together Eqs. (70) and (75), we obtain 



(to, + m) Mb (2GeV) = 102(3)(16) MeV (76) 

for the sum of the strange and light quark masses. As for 
the finer lattice, using the value m s /m = 24.4(1.5) from 
chiral perturbation theory [46], we obtain 



mf b (2GeV) = 98(3)(15)MeV 



(77) 



for the strange quark mass, which agrees well with the 
quenched lattice world average [41, 47, 48]. 

Using the GMOR inspired relation of Eq. (48) and our 
data for fp we get 



a 3 x(a) = -0.00473(28) 



(78) 



X(a) = -0.0195(12) GeV 3 . (79) 

Note that using the value f x = 0.123 GeV [49] would 
instead give x (a) = -0.0117(1) GeV 3 . 
Using ($iP) m = Zjpx(a), we finally get 

(^} m {2 GeV) = -0.0292(18)(98)GeV 3 

= -[308(6) (34) MeV] 3 . (80) 

Our results for the direct calculation of — a 3 x(m) using 
-a 3 x(m) = amJ2((P(x)P C (0)) ~ (S(x)S c (0))) (81) 

X 

are shown in Fig. 26. The quadratic fit discussed in Sec- 
tion II G, 

-a 3 X (m) = -a 3 x + B{am) + C{amf , (82) 

gives the result 

-a 3 X = 0.00448(22) , B = 0.798(4) , C = -0.14(2) . 

(83) 

The value we obtain in this way for a 3 x is compatible 
with the value in Eq. (78). The renormalized value is 

$ip) m (2GeV) = -0.0277(15)(93)GeV 3 

= -[302 (5) (34) MeV] 3 . (84) 
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Eq. (82) in the interval am < 0.132. 
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F. Direct comparison of the two lattices 

We compare in Figs. 27-29 our results for the pseu- 
doscalar and vector spectra for the finer and coarser lat- 
tices, using the Sommer scale value for the lattice spacing 
to express masses in physical units. We neglect logarith- 
mic effects in the lattice spacing and plot the mass spec- 
tra as a function of bare quark mass. It is interesting to 
observe that our results for the mass spectra on the two 
different lattices arc very similar. A qualitatively equiv- 
alent conclusion would be reached with the renormalized 
quark mass. Our results suggest that the scaling viola- 
tions for the quantities that we consider may be quite 
small. However, one has to keep in mind, that we only 
have data for two values of the lattice spacing and that 
the mass parameter p is not held fixed (for a detailed 
discussion of this point see Appendix A). 

Table I in Appendix B shows a direct comparison of 
data from the two lattices. 



IV. LIGHT BARYON SPECTRA 
A. Baryon spectra 

The increased temporal extent of our lattice with re- 
spect to that in [10] has made possible the calculation 
of the light baryon spectrum. Despite the challenges of 
limited statistics and the use of point sources and sinks, 
the lightest octet and decuplet masses were measured, 
together with those of the corresponding negative-parity 
states. For other recent determinations of the baryon 
spectrum with chiral fermions, see [20, 53-59]. 

Baryon correlation functions were constructed with the 
following interpolating operators. For the octet, we use 



0£ /2/3 = e ab M lT C l5 ^ f b 2 W c % 



(85) 



where C = 7274 is the charge conjugation matrix, a is 
a Dirac index, fi denotes flavor, and a, 6, c denote color. 
For the decuplet states, we use a notation similar to that 
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in [34] and define 

r± = (72Ti7i)/2, r = i73 . (86) 

We work in a representation of the Dirac matrices where 
74 = diag(l, 1, — 1, —1). Labeling the channels by J z , the 
decuplet operators are then given by 

D{f h = e ab M lT CT + ^)^ a=1 , (87) 
D -l/t = e abc {¥ a ^CT_^)4% , (88) 

+ (^ lT Cr + ^)^]/3, (89) 

+(^ T cr_^)^]/3. (90) 

Without loss of generality, we may assume that the 
three quark flavors are distinct and for simplicity call the 
flavors u, d, s. For the results that follow, two quarks are 
always taken to have the same mass, and assigning them 
identical flavor would merely change the normalization 
of the correlator. The two possible octet states may be 
identified with the S° and A and are given by 

E Q = O d a su + O u a sd , (91) 
A« = O d a su - Ol sd - 20 u a ds . (92) 

These give identical correlators only when all three quark 
masses arc degenerate. When propagating forward, the 
a = 1 , 2 components of these operators have positive- 
parity, while the a = 3, 4 components have negative- 
parity To extract a mass, we must project out states 
of definite parity by defining 

G ± (t) = £( |£ Q (x,t)(I^) Q ^(0,0)|0) . (93) 

X 

The mass m + of the positive-parity state is then ex- 
tracted by fitting 4e- m +* = G+(t) + G_(T - t) for 
an appropriate range of times t (where we can neglect 
the backward-propagating negative-parity contribution 
to G+(t) and vice versa). Here T is the total extent 
of the lattice in the time direction. Similarly, Ae~ m ~ t = 
G_ (t) + G + (T — t) yields the mass of the negative-parity 
state. 

For the decuplet we use 

(E*)j, = D u j ds + Df d + D d jl u (94) 

and combine correlators for the four spin states with the 
corresponding time-reversed correlators of opposite par- 
ity. The negative-parity operators are defined by replac- 
ing the Dirac index of ip* 3 in D?/* 2 * 3 according to 1 — ► 3, 
2 -» 4. 

Here we have followed convention by defining our op- 
erators in covariant form. We note, however, that one 
may also work directly with components and construct 



the spin wave functions explicitly. For example, using 
the notation |a/?7) = tabc^a^t B^ca' alternate opera- 
tors for the two spin-up octet states are 

A = (|121)-|211))/V2, (95) 
E = (|121} + |211>-2|112»/Vf$. (96) 

By construction, these have only terms with upper 
components, whereas the operators given in covariant 
form also include terms involving lower components 
(e.g. |341)). Contributions from these additional terms 
are suppressed since the lower components vanish in the 
non-relativistic limit, and the two types of operators were 
in fact found to give compatible results for both the octet 
and decuplet. For the results presented here, only the co- 
variant forms were used. 

Baryon masses were calculated at J3 = 6 with two de- 
generate quarks having each of the five lightest avail- 
able masses (ami = am 2 = 0.03,0.04,0.06,0.08,0.1) 
and for all available masses of the third quark (includ- 
ing 0.25, 0.5, and 0.75). Errors were estimated by a 
bootstrap procedure where correlators for the four or 
eight spin channels are grouped by configuration prior 
to sampling. Fitting windows were chosen based on 
plots of the effective mass M e ff = ln[C(t — a)/C(t)] for 
ami = ami = am^ = 0.03. The preferred window (used 
below for chiral fits) was determined by choosing i m , n 
such that M e ff(t m i n ) and M e ff(t m i n + a) are compatible 
to within 1 a and t max such that the error in M e ff(t max ) 
(determined by the bootstrap method) does not exceed 
30 percent. The latter criterion was made more strin- 
gent for the positive-parity decuplet state for reasons de- 
scribed below. The windows chosen on the basis of the 
lightest quark mass were then used for all other quark 
masses. 

For the octet states, the windows were 8 < t/a < 16 
for J p = | + and 6 < t/a < 8 for J p = | . Data for a 
range of windows are provided in Table X in Appendix B. 
In Fig. 30, we plot the A-like octet masses as a function 
of total quark mass. Measurements for two values of t m i n 
are shown in order to give some indication of the depen- 
dence on fitting window. Figure 31 shows the splitting 
between the two octet states as the mass of the third 
quark is increased away from ami = am2 = 0.03. 

Results for the decuplet states are shown in Fig. 32. 
The data are reproduced for a range of fitting windows 
in Table XI. We note that the masses of the heavier 
states ( J p = i , | ) exhibit a substantial dependence 
on the fitting window. In particular, the effective mass of 
the positive-parity decuplet state is seen to first plateau 
and then, at large times, continue downward toward the 
octet mass. This behavior may be explained by consid- 
ering what gives rise to the more rapid fall-off of correla- 
tors for excited states. Correlators for the octet and the 
decuplet are constructed from the same ingredients, the 
same quark propagators. The more rapid fall-off of the 
decuplet arises because of cancellations between terms. 
A fit of the correlator gives a reliable determination of 



17 



2 
1.8 
1.6 
1.4 
1.2 

i 1 
0.8 

0.6 

0.4 

0.2 










0.2 



O J P =l/2 + , t/a = 6 -16 - 



J =1/2 , t/a = 



-16 - 



J =1/2", t/a = 5 - 



A J =1/2", t/a = 6 - 8 



0.4 



0.6 



0.8 



afmj+nij+mj) 



2 
1.8 
1.6 
1.4 
1.2 

i 1 
0.8 

0.6 

0.4 

0.2 










0.2 



0.4 



O J =3/2 , t/a = 6 -16 - 
J P =3/2 + , t/a = 8 -16 - 
<■ /= 3/2", t/a = 7 -10 
A /= 3/2", t/a = 8 -10 



0.6 



0.8 



a(m ] +m 2 +m 3 ) 
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FIG. 32: Positive- and negative-parity decuplet masses for 
two fitting windows and quark masses mi = mi degenerate. 
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FIG. 33: Baryon masses with three light degenerate quarks 
were extrapolated to the chiral limit with assumed linear de- 
pendence on quark mass. 



the mass only insofar as these cancellations are not over- 
whelmed by fluctuations. At large times, the remaining 
fall-off is due primarily to the constituent masses of the 
quarks rather than the baryon mass itself. For this rea- 
son, we have constrained t max to the value used for the 
octet, giving the fitting window 8 < t/a < 16. The uncer- 
tainty associated with this choice of window is relatively 
large and of the same order as the statistical error. The 
window for the J p = | state is 8 < t/a < 10. We also 
note that the windows for both negative-parity states are 
rather small due to the early onset of fluctuations. For 
all of these reasons, our results for the negative-parity 
octet and decuplet states should be considered to have 
indicative value only. 

Finally, in Fig. 33 we plot the masses of the positive- 
parity states for light degenerate quark masses ami = 
am 2 = am 3 = 0.03,0.04,0.06,0.08,0.1. A linear extrap- 
olation to the chiral limit gives oM 8 = 0.559(24) and 
aMio = 0.690(32). In this limit, using the value of aM p 



from Eq. (18), we find 

M s /M p = 1.37(8), (97) 
Mw/M 8 = 1.23(8) . (98) 

The limited statistics of our data do not warrant a more 
complicated fitting form. Experimentally, M^/M p = 
1.21 and M A /M N = 1.31 [39]. 



B. Baryon scaling 

To investigate scaling, masses of the positive-parity 
states were also calculated at (3 = 5.85 on the 14 3 x 48 
lattice. In light of the small fitting windows that had 
been required on the finer lattice, fits for the correspond- 
ing negative-parity masses were not attempted. Here we 
find aM 8 = 0.741(28) and aM w = 1.033(55) in the chiral 
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FIG. 34: Baryon masses with light degenerate quarks at two 
values of /3, with lattice spacing set by the Sommer scale. 



limit, yielding 



Ms/M p 
M 10 /M 8 



1.23(8) , 
1.39(9) . 



(99) 
(100) 



As in the meson analysis, we make use of the Som- 
mer scale, which gives a -1 = 2.12 GeV at /3 = 6 and 
a- 1 = 1.61 GeV at = 5.85. Figure 34 shows the 
J p = i + and J p — | + states for light degenerate quarks 
on both lattices. Bare masses are rescaled by the corre- 
sponding values of a -1 . We see that the octet spectrum 
exhibits good scaling while the decuplet shows some indi- 
cation of scaling violation. Comments in Section III F re- 
garding renormalized quark masses, however, apply here 
as well. Also, as we have already observed, the decuplet 
masses suffer from uncertainty in the choice of fitting 
window and so the apparent lack of scaling should not 
be considered to have much significance. 



V. DIQUARK CORRELATIONS 
A. Diquark correlators 

In the past two years, experiments have produced indi- 
cations of bound quark systems beyond the usual quark- 
antiquark mesons and three quark baryons [60-64]. The 
most prominent example is the five particle O + (1540) 
pentaquark. The reality of such states is in question, 
since a number of attempts at verifying the pentaquark 
observations have failed [65-72]. 

Nevertheless, lattice calculations should provide a 
definitive check on whether QCD predicts such states. 
In certain models [73] of the O + (1540), which consists of 
the bound valence quarks uudds, the four quarks bind 
into two pairs of diquarks, and the two diquarks bind 
with the remaining s. 

Information on possible diquark states can be obtained 
by measuring diquark correlations on the lattice [74]. A 




FIG. 35: 3 positive-parity diquark correlation function; 
ami = 0.03, am,2 = 0.03. 



study of correlations inside baryons, using a method simi- 
lar to the one we used to investigate the qq wave functions 
in Section II D. is in progress. 

The fact that our propagators were calculated in Lan- 
dau gauge allows us also to measure quark-quark corre- 
lations directly and to fit their decay in Euclidean time 
in terms of an effective "diquark mass" [75]. Of course, 
because of the gauge fixing, one should not consider such 
a mass parameter to be the mass of a physical state. 
Nevertheless, it can produce an indication of the rela- 
tive strength of quark bindings inside diquark states. We 
consider correlations for the diquark operators 



(101) 



and 



a) (102) 



which are a 3 and 6 of color, respectively. Using these 
operators, we form four types of diquark states: (i) color 
3, spin-0, flavor 3, (ii) color 3, spin-1, flavor 6, (iii) color 
6, spin-0, flavor 6, and (iv) color 6, spin-1, flavor 3. 

Diquark correlation functions for ami = 0.03 and 
arri2 = 0.03, the lightest quark mass combination, arc 
displayed in Figs. 35-38. Since 74 was diagonal in the 7- 
matrix basis we used, upper components were combined 
with time-reversed lower components to form positive- 
parity states, while mixed upper and lower components 
were combined to form negative-parity states. Figure 39 
shows a plot of the quark correlator for input quark mass 
0.03. 



B. Diquark spectra 

Figure 40 shows the (positive-parity) diquark spectrum 
and the constituent quark masses as a function of input 
quark mass. The fitting window used was 5 < t/a < 15. 
That figure also includes a plot of twice the constituent 
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FIG. 36: 6 positive-parity diquark correlation function; FIG. 38: 6 negative-parity diquark correlation function; 
ami = 0.03, am 2 = 0.03. ami = 0.03, am 2 = 0.03. 




FIG. 37: 3 negative-parity diquark correlation function; 
ami = 0.03, am 2 = 0.03. 
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FIG. 39: Quark correlation function for input mass am 
0.03. 



quark mass and extrapolations to zero quark mass for 
the "twice quark mass" and 3 spin-0 results. Figure 41 
shows the same for the negative-parity states, with a fit 
for the lowest energy diquark state, the 3 spin-1. 

It is interesting to observe that the 3 diquark spin-0 
extrapolation in Fig. 40 is below twice the quark mass 
extrapolation and that the 3 diquark spin-0 state is sig- 
nificantly more strongly bound than the 3 diquark spin- 
1 state. Such a result is consistent with the predictions 
of diquark models [73, 76]. However, a much more de- 
tailed analysis must be done, particularly on diquarks 
within baryon states, in order for rigorous conclusions to 
be reached. 

The diquark data is shown in Tables XIV and XV in 
Appendix B. Constituent quark masses, calculated from 
fits (5 < t/a < 15) to the quark correlators, are shown in 
Table XVI in Appendix B. 
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VI. CONCLUSIONS 



In this article, we have presented results from quenched 
lattice QCD simulations using the overlap Dirac operator 
on an 18 3 x 64 lattice at p = 6.0 and a 14 3 x 48 lattice 
at (3 = 5.85. We calculated quark propagators with a 
fixed source point and a variety of quark mass values 
for 100 configurations on each lattice, and we have used 
them to evaluate meson and baryon observables, quark 
masses, meson final-state wave functions, and diquark 
correlations. 

One important result of our work is that the calcu- 
lation of quark propagators with the overlap operator 
to a chosen numerical precision has been shown to be 
feasible using available techniques, even when dealing 
with rough background gauge configurations. Such a 
calculation requires projection of low lying eigenvectors 
of H 2 = (75-Dvk) 2 - Improved algorithms or smoothing 
techniques both may help to make the calculation less 
computationally demanding [77]. 

Beyond this, our investigation validates the good chiral 
properties of the overlap operator and suggests good scal- 
ing properties between (3 = 5.85 and j3 = 6, indicating 
that the (3 = 6 results may already be close to the contin- 
uum limit. So far as the actual values of the observables 
are concerned, our results suffer from the shortcomings 
of the quenched approximation. Nevertheless, from this 
investigation and others it is clear that it should be pos- 
sible to use the overlap operator in dynamical fermion 
simulations, at the very least with a mixed action for- 
mulation. Work in that direction is beginning. We arc 
also working to extend our diquark results to diquark 
correlations in baryons with one heavy and two light 
quarks, and we expect to soon publish our results on 
non-pcrturbative renormalization, selected weak matrix 
elements, and heavy quark observables. 
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APPENDIX A: THE OVERLAP FORMALISM 
AND SIMULATION DETAILS 

The forward and backward lattice covariant derivatives 
are denned by 

V„tp(x) = ~ [U^{x)Hx + apt) - tp(x)] , (Al) 

and 



where 



75 = 75 1 



D 



(A9) 



which satisfies jI = 75 and 75 = 1. Invariance of the 
action under non-singlet chiral transformations, defined 
including a flavor group generator in Eq. (A8), forbids 
mixing among operators of different chirality [81]. There- 
fore, 



V;V(z) - - \4>{x) - Ul(x - afrUx - a£)] , (A2) 

where U )1 (x) are the gauge link fields on the lattice. Let 
Dw denote the Wilson-Dirac operator 



d w = ^7m(v m + v;)- ^v;v,, 



(A3) 
1 iii our 



with < r < 1 and < p < 2r (we used r = 
calculations) and where a is the lattice spacing. 

Neubcrgcr's overlap Dirac operator is then defined 
as [78, 79] 



D = £(1 + V)= P - (l + l5 H^l= 
a a V VHW 



where 



H = l5 {D 



w 



a 



(A4) 



(A5) 



and p is a parameter that affects the radius of rescaling of 
eigenvalues of the overlap operator in the complex plane 
as compared to the Wilson-Dirac operator. 

The combination of the Wilson gauge action and the 
(massive) overlap fermionic action is 



S 



-y 

9l V 



+ip 



1 _ I Tr (Up + U P ) 



1 — — m 
2p 



D 



*l> , (A6) 



where Up is the Wilson plaquette, go = x/6//3 is the bare 
coupling constant, the fermion fields ip an d tp carry im- 
plicit color, spin, and flavor indices, and to is a diagonal 
matrix of bare masses (tni, rn,2, ■ ■ ■ ) in flavor space. 
The Ginsparg- Wilson relation [80] 



• no additive quark mass renormalization is required, 
and the quark mass which enters the vector and 
axial Ward identities is the bare parameter to (a); 

• masses and matrix elements are affected only by 
0(a 2 ) errors and no fine-tuned parameters are re- 
quired to remove 0(a) effects; 

• and the chiral condensate (see Section II G) does 
not require subtractions of power divergent terms 
(in the chiral limit). 

The non-singlet "local" (source and sink at the same 
point x) bilinear operators we use are defined by 



O r (x)=Mx)T 



(1 



2p 



(A10) 



where Or € { S, P, Vp , A^ } correspond to T e 
{1,75,7^,7^75}- (We also use non-local "extended" op- 
erators; see Section II D.) The bilinear operators are sub- 
ject to multiplicative renormalization only, i.e. the corre- 
sponding rcnormalizcd operators are 



Or(x,fx) = lim Zr(ap) Or(x, a) 

a— >0 



(AH) 



where Zr(ap) are the appropriate renormalization con- 
stants. Since S, P and V^^A^ belong to the same chiral 
multiplcts, Zs = Zp and Zy = Za- Also, flavor symme- 
try imposes Zg = 1/Z m , with Z m defined in Eq. (37). 

Using such operators, the following two-point correla- 
tion functions (or "correlators") can be formed: 



G SS (t) = ]T(S(f,i)S c (0,0)) 



(A12) 



Gpp(t) =^<P(x,i)P c (0,0)) , (A13) 



75-D + D75 = - D l5 D 
P 



(A7) 



is satisfied by the fermionic operator of the overlap ac- 
tion, implying an exact continuous symmetry of the ac- 
tion in the masslcss limit [9]. The symmetry can be in- 
terpreted as a lattice form of chiral invariance at finite 
cutoff, 



Sip = 75^, Sip = 07 5 , 



(A8) 



Gv A p(t) = ]T<Vo.4o(z, t)P c (6, 0)) , (A14) 



G vv (t) = Y^{Vi(x,t)V^(0,0)) 



(A15) 



where we assume the quarks to be of different flavor, c 
denotes flavor conjugation, i.e. interchange of the two 
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flavors, and Vo is the symmetric lattice derivative in the 
time direction. When appropriate, such correlators are 
symmetrized around t = T/2, where T is the extent of 
the lattice in the time direction. 

Our simulations studied quenched QCD with the Wil- 
son gauge action and with Neuberger's overlap Dirac op- 
erator for lattice fermions on two different lattices. For 
both simulations, we used samples of 100 gauge config- 
urations generated by a 6-hit Metropolis algorithm with 
acceptance ~ 0.5. The configurations were separated by 
10,000 upgrades, after an initial set of 11,000 upgrades 
for equilibration. 

The first simulation was done on an 18 3 x 64 lattice, 
with (3 = 6, p = 1.4, and lattice spacing a" 1 ~ 2.0 GcV. 
We calculated overlap quark propagators for a single 
point source, for all 12 color-spin combinations and quark 
masses am = 0.03,0.04,0.06,0.08,0.1,0.25,0.5 and 0.75, 
using a multi-mass solver. 

The second simulation was done on a 14 3 x 48 lat- 
tice, with (3 = 5.85 and p = 1.6. The second, coarser 
lattice was chosen to have roughly the same volume as 
the finer 18 3 x 64 lattice, with lattice spacing a -1 ~ 
1.5 GeV, allowing for an investigation of scaling effects. 
For the 14 3 x 48 lattice, we calculated overlap quark 
propagators again for a single point source and all 12 
color-spin combinations, and with bare quark masses 
am = 0.03,0.04,0.053,0.08,0.106,0.132,0.33,0.66, and 
0.99. The largest eight of the nine quark masses cor- 
respond approximately to the eight quark masses on the 
finer lattice. 

One question that needs to be addressed when going 
to a different coupling is how to set the negative mass 
parameter p of the overlap operator. A possible and valid 
strategy is, of course, to keep the value of p constant. One 
has to keep in mind, however, that p = 1.4 was chosen 
to maximize the locality of the overlap operator at one 
specific coupling ((3 = 6.0) [27] and, from a technical 
point of view, is not optimal when going to lower values 
of 13. 

We are not restricted, however, to keep p fixed. In fact, 
since a change in p corresponds to an C(a 2 ) redefinition 
of the overlap operator (as long as we stay within the 
1-fcrmion sector of the theory), varying p{a) = po + /(a), 
where /(a) is a smooth monotonic function in a with 
lim a ^o / (a) = 0, does not change the continuum limit 
or 0(a 2 ) scaling violations of the theory. It is very rea- 
sonable to assume that choosing p(a) by demanding op- 
timum locality of the resulting overlap operator falls into 
this class of variations, and we therefore followed this 
strategy. Apart from better locality, that choice has an 
additional benefit: the resulting Hcrmitian Wilson oper- 
ator is better conditioned and therefore the number of 
eigenmodes that need to be treated exactly is smaller, 
which is particularly important at the relatively large 
physical volume we are working at. 

We used two algorithms for numerical implementation 
of the overlap operator. The first was a Zolotarev optimal 
rational function [82-84] approximation with 12 poles, as 



detailed in [11]. The second algorithm was a Chcbyshcv 
polynomial approximation [82, 83, 85, 86]. In both cases, 
we performed a Ritz projection [87] of a certain number 
ni of the lowest eigenvectors of H 2 = (j 5 Dw) 2 , whose 
contribution to the propagators was calculated directly. 
We used the Zolotarev approach for calculating 55 of the 
quark propagators for the 18 3 x 64 lattice, but numerical 
experimentation showed the Chebyshev approach to be 
about 20% faster than the rational function approach. 
We used the Chebyshev approach for the balance of the 
propagators on the 18 3 x 64 lattice and for all 100 quark 
propagators on the 14 3 x 48 lattice. 

The maximum degree of the Chebyshev polynomials 
was chosen so as to achieve a required numerical precision 
in the range of eigenvalues left over after the projection. 
We emphasize that our method of implementing the over- 
lap is exact, up to the chosen numerical precision e\ in 
the calculation of 1/ \/{H 2 ) and the maximum residue e 2 
in the calculation of the propagators. No other approxi- 
mations are involved in the calculation. We used as con- 
vergence criteria e x = 1.0 x 10 -8 and e 2 = 1-0 x 10~ 7 . For 
the 18 3 x 64 lattice, we found a projection of ni = 12 low 
eigenvectors to be adequate, leading to a maximum de- 
gree 100 ~ 500 in the expansion of the inverse square root 
into Chebyshev polynomials. However, for the 14 3 x 48 
lattice we found that, most likely on account of the in- 
creased disorder of the gauge background, the spectrum 
of H 2 generally contained many more low lying eigenval- 
ues. For some configurations, using ni = 12, as for the 
larger lattice, led to a maximum degree of the Chebyshev 
polynomials as high as ~ 4000, and in some cases led to 
loss of convergence. Therefore, for the 14 3 x 48 lattice 
we increased the number of projected low eigenvectors 
to ni = 40. This brought the maximum degree of the 
Chebyshev polynomials back into the range 100 ~ 500. 

The computations were performed with shared mem- 
ory Fortran 90 code, optimized and run on 8, 16, and 32 
processor IBM-p690 nodes at BU and NCSA. 

APPENDIX B: TABLES 

Here we present a comparison between data for (3 = 6 
and (3 = 5.85, as well as results for meson, baryon, and 
diquark spectra. 
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TABLE I: Comparison of data for the two lattices. 



On ant itv 

Vc^ Lid ±1 L 1 Lr y 




1 8 3 X 64 lattice 




1 4 3 x 48 lattice 


7 A 




L.ooooy^ t j 




1 AA"\A( 1 R\ 


A (npcptipraf p t 




29('M 




1 7(A) 


S (non-degenerate) 




0.18(8) 




0.22(4) 


a" 1 (Sommer) [24, 25] 




2.12 GeV 




1.61 GeV 


a" 1 (physical planes) 




2.19(6) GeV 




1.44(4) GeV 


a' 1 (M p ) 




1.90(4) GeV 




1.28(6) GeV 


fr/fir 




1.13(4) 




1.09(4) 


JK* 1 Jp 




l.Uo^D ) 






M K */M P 




1.09(5) 




1.07(6) 


Zf s {2GeV) 




1.44(2)(3) 




1.48(3)(16) 


(m„ + m) W (2 GeV) 




105(3)(4) MeV 




102(3)(16) MeV 


m?S(2 GeV) 




101(3)(4) MeV 




98(3)(15) MeV 


(^V) MS (2GeV) 




-[262(5) (9) MeV] 3 




-[302(5) (34) MeV] 3 


r (m s +m) MS (2GeV) 




0.265(7) (5) 




0.258(7) (3) 


r 3 W) MS (2GeV) 




-[0.684(16)(4)] 3 




-[0.78(2)(3)] 3 


TABLE II: 18 3 x 64 point sink results for PP correlator 


train 


tmax 


am 


aM 


a 3 Z/(2M) 


12 


32 


0.030 


0.2186(26) 


0.00766(64) 


12 


32 


0.040 


0.2468(22) 


0.00654(48) 


12 


32 


0.060 


0.2960(17) 


0.00560(32) 


12 


32 


0.080 


0.3396(14) 


0.00531(25) 


12 


32 


0.100 


0.3795(12) 


0.00526(21) 


12 


32 


0.250 


0.6281(9) 


0.00664(16) 


12 


32 


0.500 


0.9805(7) 


0.01094(21) 


12 


32 


0.750 


1.3046(9) 


0.01755(34) 




TABLE III: 18 3 x 64 point sink results for PP - 


5*5 correlator 




tmin 


tmax 


am 


aM 


a 3 Z/(2M) 


12 


32 


0.030 


0.2192(30) 


0.00770(70) 


12 


32 


0.040 


0.2474(24) 


0.00660(48) 


12 


32 


0.060 


0.2967(19) 


0.00567(32) 


12 


32 


0.080 


0.3403(16) 


0.00539(27) 


12 


32 


0.100 


0.3803(14) 


0.00536(24) 


12 


32 


0.250 


0.6304(10) 


0.00710(18) 


12 


32 


0.500 


0.9838(8) 


0.01206(22) 


12 


32 


0.750 


1.3084(10) 


0.01971(37) 


TABLE IV: 18 3 x 64 point sink results for VV correlator 


train 


tmax 


am 


aM 


a 3 Z/(2M) 


8 


32 


0.030 


0.511(21) 


0.00243(43) 


<s 


32 


0.040 


0.500(16) 


0.00207(30) 


8 


32 


0.060 


0.500(10) 


0.00183(18) 


8 


32 


0.080 


0.515(7) 


0.00184(13) 


8 


32 


0.100 


0.536(6) 


0.00195(11) 


8 


32 


0.250 


0.728(2) 


0.00340(10) 


8 


32 


0.500 


1.055(1) 


0.00674(14) 


8 


32 


0.750 


1.381(1) 


0.01204(24) 
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TABLE V: 18 3 x 64 extended sink results for VV correlator 



train 


traax 


am 


aM 


4 


32 


0.030 


0.446(13) 


4 


32 


0.040 


0.453(10) 


4 


32 


0.060 


0.472(8) 


4 


32 


0.080 


0.495(6) 


4 


32 


0.100 


0.519(5) 


4 


32 


0.250 


0.721(2) 


4 


32 


0.500 


1.053(1) 


4 


32 


0.750 


1.380(1) 



TABLE VI: 14 3 x 48 point sink results for PP correlator 



train 


traax 


am 


aM 


a 3 Z/(2M) 


10 


24 


0.030 


0.2513(29) 


0.0251(13) 


10 


24 


0.040 


0.2858(25) 


0.0211(10) 


10 


24 


0.053 


0.3252(22) 


0.0182(8) 


10 


24 


0.080 


0.3959(21) 


0.0157(6) 


10 


24 


0.106 


0.4557(20) 


0.0150(6) 


10 


24 


0.132 


0.5101(19) 


0.0149(5) 


10 


24 


0.330 


0.8502(13) 


0.0190(5) 


10 


24 


0.660 


1.3268(11) 


0.0357(9) 


10 


24 


0.990 


1.6488(43) 


0.0333(17) 



TABLE VII: 14 3 x 48 point sink results for PP - SS correlator 



train 


traax 


am 


aM 


a 3 Z/(2M) 


10 


24 


0.030 


0.2458(51) 


0.0232(20) 


10 


24 


0.040 


0.2821(43) 


0.0202(15) 


10 


24 


0.053 


0.3228(36) 


0.0177(12) 


10 


24 


0.080 


0.3950(29) 


0.0155(8) 


10 


24 


0.106 


0.4555(25) 


0.0149(7) 


10 


24 


0.132 


0.5106(23) 


0.0149(6) 


10 


24 


0.330 


0.8532(13) 


0.0200(6) 


10 


24 


0.660 


1.3317(11) 


0.0395(11) 


10 


24 


0.990 


1.6595(41) 


0.0420(21) 



TABLE VIII: 14 3 x 48 point sink results for VV correlator 



train 


traax 


am 


aM 


a 3 Z/(2M) 


8 


24 


0.030 


0.775(81) 


0.0151(40) 


8 


24 


0.040 


0.732(54) 


0.0104(22) 


8 


24 


0.053 


0.701(36) 


0.0076(14) 


8 


24 


0.080 


0.692(21) 


0.0060(8) 


8 


24 


0.106 


0.708(15) 


0.0058(6) 


<s 


24 


0.132 


0.732(11) 


0.0058(5) 


<s 


24 


0.330 


0.984(3) 


0.0096(3) 


8 


24 


0.660 


1.433(2) 


0.0212(5) 


8 


24 


0.990 


1.857(4) 


0.0463(20) 
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TABLE IX: 14 3 x 48 extended sink results for VV correlator 



''rain 




tmax 


am 




aivi 


A 

4 




24 




a aoa 
0.U3U 




a ^?KA/'o^^ 
O.DOU(2o) 


4 




24 




0.040 




0.645(21) 


4 




24 




0.053 




0.646(17) 


4 




24 




0.080 




0.664(12) 


4 




24 




0.106 




0.689(9) 


4 




24 




0.132 




0.719(8) 


4 




24 




0.330 




0.979(3) 


4 




24 




0.660 




1.433(2) 


4 




24 




0.990 




1.877(4) 




TABLE X: Baryon octet spectrum at f3 


= 6.0 for multiple fitting 


; windows 




am 








aM 








6-16 




8-14 


8-16 


8-18 


10-16 


0.030 


0.650(19) 




0.625(22) 


0.627(23) 


0.632(24) 


0.626(32) 


0.040 


0.680(14) 




0.658(16) 


0.658(17) 


0.662(17) 


0.653(22) 


0.060 


0.7353(97) 




0.717(11) 


0.714(11) 


0.714(11) 


0.705(14) 


0.080 


0.7854(77) 




0.7691(87) 


0.7634(88) 


0.7606(88) 


0.753(11) 


0.100 


0.8327(65) 




0.8172(74) 


0.8103(74) 


0.8052(74) 


0.7991(89) 




TABLE XI: Baryon decuplet spectrum at /3 = 6.0 for multiple fitti 


ng windows 




am 








aM 








6-16 




8-14 


8-16 


8-18 


10- 16 


0.030 


0.857(26) 




0.785(31) 


0.751(32) 


0.728(34) 


0.643(36) 


0.040 


0.863(21) 




0.805(23) 


0.775(23) 


0.757(24) 


0.694(25) 


0.060 


0.889(16) 




0.846(15) 


0.824(15) 


0.812(16) 


0.774(17) 


0.080 


0.920(12) 




0.885(12) 


0.868(12) 


0.859(12) 


0.832(13) 


0.100 


0.953(10) 




0.9222(94) 


0.9085(96) 


0.901(10) 


0.881(11) 




TABLE XII 


: Baryon octet spectrum at (3 


= 5.85 for multiple fitting windows 




am 








aM 








6-11 




7-10 


7-11 


7-12 


8-11 


0.030 


0.805(35) 




0.796(43) 


0.813(48) 


0.817(52) 


0.839(70) 


0.040 


0.840(23) 




0.833(27) 


0.840(30) 


0.838(33) 


0.849(41) 


0.053 


0.879(15) 




0.872(17) 


0.873(19) 


0.868(21) 


0.873(25) 


0.080 


0.9526(94) 




0.946(11) 


0.941(11) 


0.935(12) 


0.935(13) 


0.106 


1.0196(78) 




1.0130(93) 


1.0065(88) 


0.9996(88) 


0.998(10) 


0.132 


1.0842(69) 




1.0772(82) 


1.0700(76) 


1.0630(74) 


1.0602(84) 




TABLE XIII: 


Baryon decuplet spectrum at /3 — 5.85 for multiple fitting windows 




am 








aM 








6-9 




7-9 


7-10 


7-11 


8-10 


0.030 


1.131(52) 




1.106(77) 


1.123(90) 


1.131(98) 


1.17(17) 


0.040 


1.140(38) 




1.110(53) 


1.118(60) 


1.123(65) 


1.135(94) 


0.053 


1.153(28) 




1.121(37) 


1.121(40) 


1.121(41) 


1.118(54) 


0.080 


1.190(18) 




1.157(22) 


1.148(22) 


1.142(21) 


1.133(26) 


0.106 


1.233(13) 




1.201(15) 


1.189(15) 


1.180(14) 


1.171(16) 


0.132 


1.279(11) 




1.248(12) 


1.235(11) 


1.224(10) 


1.217(12) 
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TABLE XIV: Degenerate (mi = 7712) diquark spectrum, positive-parity. 



color 


spin 


am 


aM 


3 





0.030 


0.430(13) 


3 





0.040 


0.453(11) 


3 





0.060 


0.494(8) 


3 





0.080 


0.531(7) 


3 


1 


0.030 


0.551(10) 


3 


1 


0.040 


0.560(8) 


3 


1 


0.060 


0.585(6) 


3 


1 


0.080 


0.612(5) 


6 





0.030 


0.548(16) 


6 





0.040 


0.555(14) 


6 





0.060 


0.579(11) 


6 





0.080 


0.609(10) 


6 


1 


0.030 


0.542(15) 


6 


1 


0.040 


0.552(12) 


6 


1 


0.060 


0.582(9) 


6 


1 


0.080 


0.616(8) 



TABLE XV: Degenerate (mi = m?) diquark spectrum, negative-parity. 



color 


spin 


am 


aM 


3 





0.030 


0.795(61) 


3 





0.040 


0.799(53) 


3 





0.060 


0.818(43) 


3 





0.080 


0.842(37) 


3 


1 


0.030 


0.551(29) 


3 


1 


0.040 


0.560(23) 


3 


1 


0.060 


0.585(17) 


3 


1 


0.080 


0.612(14) 


6 





0.030 


0.898(104) 


6 





0.040 


0.918(88) 


6 





0.060 


0.935(67) 


6 





0.080 


0.947(55) 


6 


1 


0.030 


0.542(38) 


6 


1 


0.040 


0.552(31) 


6 


1 


0.060 


0.582(23) 


6 


1 


0.080 


0.616(20) 



TABLE XVI: Constituent quark mass 



am 


aM 


0.030 


0.229(5) 


0.040 


0.235(5) 


0.060 


0.248(4) 


0.080 


0.261(4) 
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